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ABSTRACT 

We have used a many-body model Hamiltonian to study the nature of the magnetic ground 
state of hetero-binuclear complexes involving rare-earth and copper ions. We have taken 
into account all diagonal repulsions involving the rare-earth 4f and 5d orbitals and the 
copper 3d orbital. Besides, we have included direct exchange interaction, crystal field 
splitting of the rare-earth atomic levels and spin-orbit interaction in the 4f orbitals. We 
have identified the inter-orbital Af repulsion, U// and crystal field parameter, Ay as the 
key parameters involved in controlling the type of exchange interaction between the rare 
earth Af and copper 3d spins. We have explored the nature of the ground state in the 
parameter space of U/y, Af, spin-orbit interaction strength X and the Af filling ny. We 
find that these systems show low-spin or high-spin ground state depending on the filling of 
the Af levels of the rare-earth ion and ground state spin is critically dependent on Uyy and 
Af. In case of half-filling (Gd(III)) we find a reentrant low-spin state as U/y is increased, 
for small values of Af, which explains the recently reported apparent anomalous anti- 
ferromagnetic behaviour of Gd(III)-radical complexes. By varying U/y we also observe 
a switch over in the ground state spin for other fillings . We have introduced a spin- 
orbit coupling scheme which goes beyond L-S or j-j coupling scheme and we find that 
spin-orbit coupling does not significantly alter the basic picture. 

1 Introduction 

Electron spin pairing and the associated electron exchange between atoms is the corner- 
stone of chemical bond formation. Hence most of the molecules have closed shell ground 
state [IJ and the synthesis of high spin molecules as well as molecular magnets has been a 
challenge. The simplest idea to synthesize a high spin molecule is to start with a high spin 
carrier. In this regard lanthanides, and in particular, Gd 3+ which has a Af 1 configuration 
is an obvious choice. Mononuclear complexes of Gd 3+ are already known. Systems with 
ground state spin higher than | can be built up with multinuclear complexes involving 
Gd 3+ and other rare-earth or transition metal ions. Studies of such heteropolymetallic 
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compounds have also been of great interest in modeling metaloenzymes and in under- 
standing their magnetic properties [Q]. Indeed, the interaction between two nonequivalent 
metal centers can lead to some unique features which are not encountered in the com- 
plexes with same kind of spin carriers. 

In the last few years, there have been several studies of heteropolymetallic molecu- 
lar coordination complexes of lanthanides and Cu(II) or other ions because of their lu- 
minescence and interesting magnetic properties [J3p. Some of these complexes exhibit 
a high spin ground state corresponding to a ferromagnetic alignment of magnetic mo- 
ments of the individual ions, a fact which is important for designing building blocks for 
molecular magnets. Gatteschi and coworkers found that in case of Gd(III)-Cu(II) 
binuclear and Cu(II)-Gd(III)-Cu(II) trinuclear complexes, the Gd(III)-Cu(II) interaction 
is ferromagnetic in nature, independent of the structure of the complex or the ligand 
molecule. Subsequently, a systematic effort to synthesize and study binuclear systems 
between Cu(II) and other rare-earth ions was undertaken [Q-^] . These studies indicate a 
general trend for the ground state spin of the complex as we traverse the lanthanide series. 
The interaction between Cu(II)-Ln(III) is ferromagnetic for Gd(III) and other rare-earth 
ions with greater than half-filled 4/ shell such as Dy(III) and Er(III). Below half filling, 
this interaction is antiferromagnetic and leads to low-spin ground state for systems such as 
Cu(II)-Ce(III) and Cu(II)-Eu(III) [0]. The situation remains the same when the Cu(II) is 
replaced by organic ligands like nitronyl nitroxide triazole which are essentially spin- 1/2 
carriers. 

As a possible mechanism for the observed ferromagnetic alignment of the spin of the 
rare-earth ion with the spin of the copper ion, it was pointed out [ |T0[ ] that the stabilization 
of the high spin state arises from the coupling between the ground state configuration and 
the excited state configuration arising out of the transfer of the electron from the singly 
occupied copper 3d orbital to an empty 5d orbital of the rare-earth ion. This mechanism 



is in the spirit of the superexchange mechanisms proposed by Goodenough [ J13D in 1963. 
Similar mechanism is invoked recently by Tchougreeff [|T4p to explain the ferromagnetic 
ordering in decamethylferrocenium tetracyanoethenide, by Kinoshita and coworkers [ [l"5l ] 
to explain the intermolecular ferromagnetic coupling in paranitrophenyl nitronyl nitrox- 
ide and by Girerd and coworkers [[16]] to interpret the magnetic properties of /u-oxo Mn(II) 
compounds. Indeed, the interplay of direct exchange which favours parallel spin align- 
ment on a given site (Hund's rule) and kinetic exchange which favours antiparallel align- 
ment of spins on bonded pair of sites (by virtue of Pauli principle) is what determines the 
spin state of a system [|T7|]. 

Kahn and coworkers modeled the superexchange process using an extended Hiickel 
model and obtained a value for the septet-nonet gap from the computed Jod-Cu in fairly 
good agreement with the experimental result. But in this model strong electron-electron 
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repulsions among the Af electrons arising out of the contraction of the lanthanide Af 
orbitals as well as between Af electrons and the 5d electrons have not been taken into 
account. These electron repulsions, according to the recent reports seems to play an 
important role in explaining the photoemission spectra of rare-earth compounds [JT8p . 
Another factor which was not accounted for in the above model is the splitting of the 
free ion Ln(III) ground state due to the ligand field effect, on the grounds that spin-orbit 
coupling term is an order of magnitude larger than crystal field term QT9|] . The need of 
a detailed microscopic model taking these factors into account has also been felt since 
some exceptions to ferromagnetic coupling in Gd(III) complexes were recently reported 
for Cu(II) complexes and nitronyl nitroxide radicals [^Q|-^2|]. To explain this unexpected 
antiferromagnetic coupling, suggestions were put forward that the observed coupling is 
actually the sum of two contributions [[22|], one from the direct overlap of the magnetic 
orbitals of the ligand with the / orbitals, which leads to antiferromagnetism and the other 
from the overlap with the s and d orbitals, which results in ferromagnetism. But the rare- 
earth / orbitals are highly contracted to the core which prompted Kahn et. al. [ |lTj| ] to 
neglect the transfer between the Af orbitals of Ln(III) and 3d orbitals of Cu(II). 

In this paper we present our studies of a detailed microscopic model in which we have 
taken into account explicit electron-electron correlations, crystal field effects and spin- 
orbit effects in the Lanthanide ion. We have spanned the parameter space of the model 
and have obtained phase diagrams for different fillings of the Af orbitals to elucidate the 
role played by the different parameters in controlling the ground state spin. The paper is 
organized as follows. In section 2, we introduce the many-body model Hamiltonian. We 
also briefly describe the computational schemes employed in solving the model Hamil- 
tonian. In section 3 we discuss the results obtained by our computations on the model. 
Finally, we summarize our results in section 4. 



2 Many Body Model and Computational Scheme 

The remarkable feature of these systems is that the parallel or anti-parallel alignment of 
spins in the ground state is independent of the details of structure. This led Kahn et. al. 
to propose a general mechanism by the virtue of which these systems attain either the 
high spin or the low spin state. The contraction of the Af orbitals towards the nucleus 
means that they are very weakly delocalized towards the Cu 3d orbitals and thus, there 
is a vanishingly small overlap between the Gd Af and Cu 3d orbitals. This also means 
that any interaction between the ground state configuration and the charge transfer con- 
figuration corresponding to a 3 d — Af transfer cannot stabilize the low spin 5 = 3 state. 
So, Kahn et. al. neglected these 3d — Af interactions and considered charge transfer to 
the higher 5d level, ie, the 3d electron being transfered to an empty 5d orbital of the 
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Gd(III) ion. This process results in the S = A high spin state being stabilized due to a 
strong Hund's coupling between the Af and 5d levels. This process is illustrated in Fig. 
1 . The mechanism outlined by Kahn et. al. is however, still in a non-interacting picture 
since it does not take into account explicit Coulomb repulsions among the electrons in 
different orbitals. Another factor missing in the model is the effect of the crystal field 
interaction on the various orbitals involved. In order to take these effects into account and 
understand how they modify the ground state spin in these systems, we have studied these 
systems by employing an interacting Hamiltonian. This Hamiltonian is an extension of 
the Pariser-Parr-Pople model [ ]TT] ] to more than one orbital per atom and includes direct 
exchange interaction terms between electrons in orbitals on the same site. This model has 
been greatly successful in describing ground state and low-lying excited states of various 
strongly correlated systems. 

When we have more than one orbital per site, the ZDO (zero differential overlap) 
approximation is not valid and we should include electron repulsion integrals when two 
electrons are in different orbitals on the same atom. In this case, we will have several kinds 
of nonzero two electron integrals on the same site. The interaction part of the Hamiltonian 
in this modified ZDO approximation is given by, 

H int = 1 /2 52 52 t W I Mv] (Ewt/Ejvjv' - 8 ! - ; -8 A / v £ ! - jUjV / ) ( 1 ) 

i,j ai,a/,v,v' 

Here, i, j are site indices and /j, v correspond to the orbital on a particular site. The 
two electron integrals [^y|jvJv'] contains contributions from various kinds of situations. 
The case \i = fj = v = v' and i = j corresponds to the electron-repulsion integral for 
the electrons in the same orbital of a single site. The most general case for the repulsion 
integral on the same site is realized when there are at least four orbitals per site and 
fj ^ // ^ v ^ V and i = j. For the / orbitals on the same site, we have assumed that 
(i) the intra-orbital repulsion integral [f^f^f^] is the same for all the seven / orbitals, 
independent of /j, and given by Uf, (ii) the inter / orbital repulsion integral [f^f/j\fvfv] is 
also independent of the particular pair of / orbitals /j, v and given by U ff and (iii) that 
there is also a single direct exchange integral Jff between two / orbitals corresponding to 
the integral [f^fv \fufv] ■ We have further assumed that all other electron repulsion integrals 
involving the / orbitals are zero. The direct exchange integral although smaller than the 
electron repulsion integrals in (i) and (ii) by about an order of magnitude, is crucial to the 
study of magnetic gaps. Indeed, the small value of the exchange integral is also reflected 
in the small value of the magnetic gaps. For our rare-earth system, in most cases, we have 
considered all the seven Af orbitals and one 5d orbitals for the rare-earth metal ion and a 
single 3d orbital for the Cu ion. Inclusion of 5d orbital introduces corresponding electron 
repulsion parameters U^, U fj and Jfj. Electron repulsion integrals of the form [i^lj'viv] 
which in our case corresponds to repulsion integral V/y between an electron in the Af 
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orbital of the rare-earth and an electron in the 3d orbital of the Cu ion is parametrized 
using Ohno interpolation scheme [P"2|]. All other inter-site electron repulsion integrals are 



neglected. 

Using the above parameters, the model Hamiltonian is written as , 

7 

H = £ Af f h f + A 3d h 3d + A 5d h 5d 

f=i 

i 

+ ^d-4f £ {Ef,3d + H. C.) + t 3d - 5d (E3d,5d + H. c.) 
/=! 

+ Y £ fy&f ~ ^ + ~Y hl,d ^ 3d ~ + ~Y ^ 5d ~ ^ 

7 7 

/>/'=! /=i 

7 

+ v 3d-4f 52 "/ "3d + V 3 d-5d n 3d h 5d 
/=! 

7 

+ % E [(£/,/' £/,/' + £/,/' E f,f +H.c.) - hf - h f ] 

/>/'=! 
7 

+ E [( E f,5dEf,5d + E 5dJ E f)5d + H. c.) - h f - h 5d ] (2) 

/=! 

where we define the electron hop operator Ejj which hops an electron of either spin from 
orbital/ to orbital i as = ^ aJ a 7 ; a] G (d ia ) creates (annihilates) an electron with 
spin o in the i th orbital. The first two lines in Eqn. 2 corresponds to the noninteracting 
part of the Hamiltonian. The first line corresponds to the site energy contribution. Here, 
we assume that the degeneracy of the Af levels are lifted uniformly by the nonsymmetric 
crystalline field, leading to equally spaced levels with an energy level spacing Af. The 5 d 
orbital of the rare-earth and the 3d orbital of the Cu are located at energies A$ d and A 3d 
respectively. The second line in the Hamiltonian corresponds to the transfer part between 
the individual Af orbitals and the Cu 3d orbital, with a transfer parameter t 3d -4f and a 
transfer between the rare-earth 5d orbital and the Cu 3d orbital with a transfer parameter 
hd-5d- Here, we assume that the orbitals on the same site are eigenstates of the site 
Hamiltonian and therefore do not mix. The third line in the Hamiltonian corresponds to 
Hubbard electron repulsion terms, Uf, U$d and U 3d being the repulsion parameters for 
the orbitals we have considered. The fourth line involves inter-orbital interactions on the 
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rare-earth site with U ff for repulsion between electrons in different Af orbitals and U fd 
for repulsion between electrons in the Af and 5d orbitals at the rare-earth site. The fifth 
line in the Hamiltonian corresponds to the inter-site interaction between electrons on the 
rare-earth site and those in the Cu 3d orbital. The parameters V<y_4/ and V^d-sd are 
obtained using Ohno parametrization with intersite distance of 2A. The last two lines in 
the Hamiltonian correspond to the exchange term involving orbitals on the rare-earth site. 
The exchange term has both diagonal and off-diagonal contributions, the latter combing 
from two successive electron hops. 

We set up the matrix for the above Hamiltonian in the Valence Bond (VB) basis. 
The VB basis being eigenstates of the total spin operator is ideal for this study since the 
total spin of the targeted state is known exactly. The VB method is discussed extensively 
in [|23p. The resulting Hamiltonian matrix is non-symmetric and is diagonalized using 



Rettrup's algorithm [g4|. The low energy states are determined for a range of values of 
the parameters in the Hamiltonian. These studies were carried out at various Af orbital 
fillings, corresponding to different lanthanide ions. 

We have also carried out calculations that includes spin-orbit interactions between the 
Af orbital angular momentum and the electron spin. The spin-orbit contribution is given 
by, 

Hs.o. = £ £(r<) ■ s(i) (3) 

i 

where ^(r ; ) is given by 

Ltn L c L r or 

with usual notations. Integrating over the radial part leads to the spin-orbit interaction 
parameter X and Hs.o. can be written in the second quantized form as, 

Hs.o. = I L W l ( l+1 ) ~ ™(m+l) «l +1) p«m,a 

m 

+ y/ 1(1 + 1) - m(m- 1) a] n _ l a a m # 

+ m(a t ma a m ,a ~ a m $)] (5) 

where, / takes the value 3, corresponding to the / orbitals and m varies from -3 to 3. In 
the presence of spin-orbit interactions, the total spin of the system is not conserved. Thus, 
we cannot employ the VB method for their solution. Instead we employ a constant M s 
basis, typified by Slater determinants, to set-up the Hamiltonian matrix. The Hamiltonian 



matrix in this case is symmetric and Davidson algorithm [ |25| ] gives the low-lying states 
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Since the states are now labelled by the total angular momentum J = L + S, we compute 
the expectation value of S 2 and use the nearest integer or half-integer value S such that 
the expectation value is nearest to S(S+ 1). This allows identifying the approximate total 
spin of the states. 

3 Results and Discussion 

The Hamiltonian in Eq. (2) has several parameters and exploring the phase space for 
all possible filling of the Af rare-earth orbital is very compute intensive. In order to 
meaningfully explore the parameter space, we have firstly identified the most sensitive 
parameters and fixed the remaining parameters at generally accepted values. We have 
considered three different models for the Ln-Cu system. These correspond to (a) ignoring 
the 5d levels on the rare-earth ion, (b) including all the five 5d orbitals on the rare-earth 
ion and (c) including only one 5d orbital on the rare-earth ion. In all these cases, we 
have found that the inter-orbital electron repulsion strength , Uff, in the Af orbital and 
the crystal-field splitting of the Af orbitals Ay are the parameters to which the spin of 
the ground state is most sensitive. The remaining parameters have been held fixed at the 
following values. The value of the exchange constants Jff and Jfd are fixed at O.Wf 
and 0.2Uf , respectively. The U^d on copper is fixed at 5eV, t^d-Af at O.leV and in all 
the models and in models (b) and (c) t^dsd was fixed at 0.5eV, U$d at 5eV and Ufd at A 
eV. In what follows, we discuss our results, first for the half-filled Af shell and then for 
rare-earth ions with other Af occupancies, for various values of the parameters Uff and 
Af. In the last subsection we explore the effect of spin-orbit interactions on the "spin" of 
the ground state for these parameter values. 

3.1 Af 1 system 

For the half-filled Ln ion in model (a) where the 5d orbitals on the rare-earth site is 
ignored, we always obtain a low-spin ground state, for all values of Uf and Af that we 
scanned. Changing the parameters t^d-Af and U^d within reasonable limits still leaves the 
system in a low-spin ground state. This is because, the path way for derealization of 
the electrons in the high-spin state is blocked (Fig. 2a), depriving the system of kinetic 
stabilization. Thus, this model does not produce a high-spin ground state for any value 
of Af or Uff, although experimentally most 4/ 7 -radical or Cu(II) systems are found in 
high-spin ground state. 

In the presence of the empty 5d orbitals (Fig. 2b), the Hunds rule stabilizes the state in 
which the Af and 5d electrons have the same spin alignment. This in turn implies that the 
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high-spin state in which the unpaired electron in the Cu 3d orbital has the same orientation 
as the spins in the Af orbital has lower energy than the low-spin state. Our studies on 
model (b) is therefore expected to produce a high-spin ground state. However, what is 
quite unexpected is that we find a very interesting variation in the ground state spin; the 
ground state is a high-spin state for small Af (< O.leV) and l.5eV < Uff < A.5eV at 
Uf = lOeV. For other values of Uff, the ground state is in a low-spin state. Thus, in the 
quantum phase diagram (Fig. 3), there is a re-entrant low-spin phase. For small values of 
U ff the virtual low-spin state in which the intermediate rare-earth ion configuration is Af s 
has lower energy than the configuration Af 7 5d l and this would lead to a low-spin ground 
state of the complex. At intermediate values of Uff, the virtual state Af 7 5d l has a lower 
energy than the 4/ 8 virtual state and we should expect the ground state to be a high spin 
state. For large values of Uff, the ground state configuration of the rare-earth ion is no 
longer 4/ 7 but Af 6 5d l . In this case, the hopping between the 5d and 3d orbitals favours 
antiparallel alignment of spins. The 3d to Af hop prefers parallel alignment on energy 
consideration but there is only one channel for this process as against seven channels for 
antiparallel alignment. Therefore, the ground state switches back to low-spin state for 
high values of Uff. 

We also find that the high-spin state is quite stable and continues to be the ground state, 
even when the *3</-4/ transfer integral is increased to unto twice its value of 0. leV we have 
used in all our studies. Furthermore, we find that the Af parameter which corresponds to 
the crystal field splitting of the Af levels also plays a significant part; when Ay is small, the 
cross over from high-spin ground state to low-spin ground state occurs at a higher value 
of U ff. Above a critical value of Af, the system is found only in the low-spin ground state 
(Fig. 3). 

This picture is not changed in any significant way, in going to model (c) where only 
one 5d orbital is used instead of the five 5d orbitals on the rare-earth considered in model 
(b). In view of this, in cases away from half-filled Af levels, we have considered only a 
single 5d orbital (model (c)) for reasons of computational feasibility. 

3.2 Systems with non-half-filled Af shell 

We have studied systems with various Af orbital fillings using the Hamiltonian in Eq. 2 
and the general result we find is that when Uff is small, systems with fewer than seven 
Af electrons have a high-spin ground state while those with more than seven Af electrons 
have a low-spin ground state (Fig. 6). This result follows even in model (a) where we 
have neglected the participation of the 5d orbitals in derealization. This result is quite 
simply understood from the direct exchange interaction term. In the less than half-filled 
case (Fig. 4), an electron hop from the Cu 3d orbital to the rare-earth Af orbital can result 
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in a high-spin or a low-spin rare-earth ion. If the unpaired electrons on the rare-earth have 
the same alignment as the hopping electron, we obtain the high-spin state, otherwise the 
resulting state of the rare-earth ion is a low-spin state. The direct exchange interaction 
within Af orbitals ensures that the virtual high-spin state of the ion has lower energy than 
the virtual low-spin state. Therefore, the ferromagnetic alignment of the spins on the 
Ln-Cu system is favoured over antiferromagnetic alignment in the less than half-filled 
case. In the more than half-filled case, electron-hop is only permitted when the hopping 
electron has an alignment opposite that of the rare-earth spin (Fig. 5), because of Pauli 
exclusion principle. Thus, in the more than half-filled case, we find a low-spin ground 
state. Indeed, this picture is not significantly changed when we go to model (c) where 
we have the copper 3d electron having an additional derealization path-way via the rare- 
earth 5d orbital. What we also note is that the t^d-Sd does not have a significant effect, 
unlike in the half-filled case. 

We should expect the crystal field splitting parameter, Ay, of the Af orbitals to play 
an important role in determining the ground state of the Ln-Cu system, since in case 
of large Ay, the rare earth ion ground state would not be a high-spin ground state. For a 
similar reason, we should also expect the inter-orbital repulsion, Uff also to determine the 
magnetic state of the dimer. In Fig. 6, we present some typical quantum phase diagrams 
for the less than half-filled case as well as the more than half-filled case. In the less 
than half-filled case, for large Uff, we obtain a low-spin ground state since the energy 
difference between high-spin and low-spin rare-earth virtual states obtained on electron 
transfer, is not significant. However, number of low-spin virtual states are more than 
the high-spin virtual states and we obtain a low-spin ground state for large Uff. The 
experimental system with less than half-filled Af shell are all known to be in the low-spin 
(afm) ground state. Thus, these systems all correspond to the large Uff limit. The critical 
value of Uff decreases with increase in Af occupancy from 1 to 7. 

In Fig. 6, we show some typical quantum phase diagrams for the more than half-filled 
case. Here again we note a quantum phase transition, from the low-spin ground state to 
the high-spin ground state, at large values of Uff. This transition is also qualitatively 
understood in the following way. At large Uff, the low-energy exchange path way is via 
the empty 5d orbitals. The direct exchange interaction between the 5d and Af orbitals on 
the rare-earth ensure that the virtual states corresponding to high-spin state has a lower 
energy than those corresponding to low-spin state. Thus, the dimer switches to a high- 
spin ground state with increasing Uff. All the experimental systems known in this case 
correspond to high-spin ground state. Thus, we are in the large Uff limit in real systems 
even in systems in which the rare earth ion has more than half-filled Af shell. This is also 
intuitively expected since the Af shell is chemically well isolated in rare earth compounds. 
The critical value of Uff in this case increases with increase in the Af occupancy beyond 
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half- filling. 



3.3 Role of spin-orbit interactions 

It is well known that the spin-orbit interactions in rare-earth ions are important and often 
larger than the crystal field effects, we have considered only the spin orbit interaction for 
electrons in the Af orbital of the rare-earth and have ignored it in the 5d orbital of the 
rare-earth as well as in the 3d orbital of the copper. Our method of calculation takes into 
account spin-orbit interactions at the same level as all other interactions in the Hamilto- 
nian in Eq. 2, since we set-up and solve the total Hamiltonian (Hamiltonians in Eq. 2 plus 
Eq. 3). Thus, our studies go beyond the usual L-S or j-j coupling schemes in which either 
the spin-orbit interaction is treated perturbatively within a given manifold of L and S or 
the electron repulsions are treated perturbatively over the individual spin-orbit coupled 
states. We have used values of X in the range to 1.0 eV, for the Uff and Af values [|7p. 
In the presence of spin-orbit interactions, total spin is no longer a good quantum number. 
Thus we cannot strictly classify the eigenstates by their spins as low-spin or high-spin 

Art 

states. However, we find that the expectation value of S in the low-lying states, in all the 
cases we have studied, is close enough to allowed S values in the absence of spin-orbit 
interactions. This aids in labelling the states as 'high-spin' or 'low-spin' states. 

What we find is that the role of the spin-orbit interaction in determining the ground 
state magnetic moment of the Ln-Cu dimers is negligible. In Fig. 7, we show the quantum 
phase diagram for three representative cases corresponding to less than half-filled, half- 
filled and more than half-filled Af shells, for Ay = as a function of X and Uff, with 
all other parameters being the same as in Figs. 3 and 6. We find that the Uff value 
for which the magnetic transition occurs in the ground state is hardly affected by spin- 
orbit interactions. Only the re-entrant low-spin state occurs for a slightly smaller value of 
U ff when spin-orbit interaction is sufficiently strong. We have also explored the role of 
spin-orbit interaction when the crystal field splitting of the Af orbitals Ay is varied near a 
critical value of Uff. Here we assume that the splitting of the /-orbital manifold by the 
crystal field does not quench the orbital angular momentum of the resulting state. In Fig. 
8 is shown the phase diagram as a function of X for the case of the rare-earth electron 
configuration Af 10 . We have chosen this configuration as we found the strongest effect 
due to spin-orbit interactions in this case. We find in Fig. 7 that high-spin state crosses 
over to the low- spin state at a smaller value of Ay when the spin-orbit interaction strength 
X is increased. Thus, we conclude from our studies that the magnetic state of the cluster 
is only weakly influenced by spin-orbit interactions. 
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4 Summary and Conclusions 



In this paper, we have studied the magnetic ground state of the hetero bi-metallic organic 
molecular magnets, with a rare-earth ion and a copper ion. These systems have been 
modeled with a many-body Hamiltonian, which includes the direct exchange interaction 
and the repulsion between the orbitals on the same site. The systems are shown to exhibit 
a low-spin or a high- spin ground state depending on the filling of electrons in the rare 
earth Af levels, and also critically on the inter-orbital Af electron repulsions, Uff and 
the splitting of the degenerate Af orbitals in the presence of the crystal field, Af. We 
have obtained a "phase diagram" in the space of these parameters. We find at half-filling 
a re-entrant low- spin state as Uff is increased for small values of Af. There is also a 
transition in the ground state spin at other fillings, again for small values of Af, when 
Uff is increased. The general phase diagram of the ground state spin is shown in Fig. 
9, as a function of filling and Uff. This picture is not altered by the inclusion of spin- 
orbit interactions. To obtain the experimentally observed interaction, corresponding to 
low-spin ground state below half-filling and high-spin ground state above half-filling of 
the Af shell, we need to assume fairly large inter-orbital repulsions in the Af shell. Our 
studies show that in systems where the parameters are near the phase boundaries, small 
changes brought about by the modification of the ligands can lead to drastic changes in 
the magnetic properties of the complex. 
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Figure 1 : A schematic for the mechanism which leads to the stabilization of the high spin 
state in the system. Ln 3+ is the rare-earth (Lanthanide) ion. Transfers between the Cu 2+ 
3d and the Ln 4f and 5d levels are indicated by the arrows. The 3d— >4f transfer leads to 
a higher energy virtual state than 3d— >5d transfer. 
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Figure 2: A schematic diagram showing the various possible virtual states obtained by 
electron transfer from 3d orbital of Cu(II) to orbitals of Ln(III), when Af is half-filled - 
(a) without 5d levels, (b) with all the 5d levels taken into consideration. 
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Figure 3: Phase diagram in the space of parameters U/f and Ay, for the case of rare-earth 
4f levels being half-filled. Other parameter values are as follows: Uf =10 ev, t?, d -4f=0.1 
ev, t 3d - 5d =0.5 ev, U 5d =5 ev, U 3d =5 ev, U fd =4 ev, J/f=0.1 ev and J fd =0.2 ev 
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Figure 4: A schematic diagram showing the various possible virtual states obtained by 
electron transfer from 3d orbital of Cu(II) to orbitals of Ln(III), when 4f is less than 
half-filled. 
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Figure 5: A schematic diagram showing the various possible virtual states obtained by 
electron transfer from 3d orbital of Cu(II) to orbitals of Ln(III), when 4/ is more than 
half-filled. 
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Figure 6: Phase diagram in the space of parameters U/f and Af, for the case of rare-earth 
4f levels being away from half-filling. Other parameter values are same as in Fig. 3 
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Figure 7: Phase diagram in the space of parameters Uff and A4/, for the case of rare-earth 
4f levels at less than half-filling, half-filling and more than half-filling. Other parameter 
values are same as in Fig. 3 
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Figure 8: Phase diagram in the space of parameters X4/ and Ay for the case of rare-earth 
Af levels being away from half-filling. Other parameter values are same as in Fig. 3 
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Figure 9: A phase diagram for the Ln(III)-Cu(II) system, between parameters U/f and 
tif, nf being the 4f level filling. Ay value was fixed at 0.1 ev. Other parameter values are 
same as in Fig. 3 
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